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Differential rotation and meridional flow in the solar 
supergranulation layer: Measuring the eddy viscosity 
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Abstract. We measure the eddy viscosity in the outermost layers of the solar convection zone by comparing the rotation law 
computed with the Reynolds stress resulting from f-plane simulations of the angular momentum transport in rotating convection 
with the observed differential rotation pattern. The simulations lead to a negative vertical and a positive horizontal angular 
momentum transport. The consequence is a subrotation of the outermost layers, as it is indeed indicated both by helioseismology 
and the observed rotation rates of sunspots. In order to reproduce the observed gradient of the rotation rate a value of about 
1.5 x 10 13 cm 2 /s for the eddy viscosity is necessary. Comparison with the magnetic eddy diffusivity derived from the sunspot 
decay yields a surprisingly large magnetic Prandtl number of 150 for the supergranulation layer. The negative gradient of the 
rotation rate also drives a surface meridional flow towards the poles, in agreement with the results from Doppler measurements. 
The successful reproduction of the abnormally positive horizontal cross correlation (on the northern hemisphere) observed for 
bipolar groups then provides an independent test for the resulting eddy viscosity. 
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1. Introduction and observations 

Over the last years, while observing the solar oscillations on 
longer timescales and with higher precision, it has become ev- 
ident that these oscillations play an important role in under- 
standing the solar interior, bearing more or less the only in- 
formation from the deeper parts of the sun, which cannot be 
probed otherwise. Helioseismology reveals a maximum of the 
angular velocity at all latitudes rather close to the surface, as 
shown in Fig. |(Howe et al. 2000). 

It is known, on the other hand, that sunspots rotate faster 
than the solar surface plasma by about 4 % or 80 m/s at 
all latitudes. Such a clearly verified subrotation of the outer- 
most layer of the convection zone is easiest understood as a 
result of angular momentum conservation of fluid elements 
with purely radial motions. But in this domain of the solar 
convection zone the velocity field is dominated by horizontal 
motions. Fluctuating fields with predominantly horizontal in- 
tensity should produce superrotation rather than subrotation. 
There is, however, another strong argument for considering the 
exceptional behavior of the horizontal motions in more detail. 

Ward (1965) was the first to consider the horizontal cross- 
correlation of the proper motions of sunspot groups, the faster 
of which tend to move toward the equator. He found 



on the northern hemisphere. More recent observations found 
smaller, but always positive values (Gilman & Howard 1984; 
Nesme-Ribes et al. 1993; Komm et al. 1994, see an overview by 
Meunier et al. 1997). This result has a strong implication for 
theory confirming the existence of the positive H coefficient in 
the expression d3J) for the horizontal Reynolds stress. We shall 



CD 

o 

o 
DC 



440 
420 



380 
360 





0° 




— " 15°^ 




" 30° \ 




^w-"""" 45° i 




^ - — „N 

60 







0.1 (deg/day) 2 « 2 x 10 7 cm 2 /s 2 



(1) 



Send offprint requests to: G. Rudiger 



0.50 0.60 0.70 0.80 0.90 1.00 
r/R 

Fig. 1. The internal rotation of the Sun as found by helioseis- 
mology. Image: NSF's National Solar Observatory 

demonstrate the close relation between the negative radial gra- 
dient of the rotation rate and the positive sign of the horizon- 
tal cross correlation by solving the Reynolds equation on the 
basis of new data from hydrodynamic simulations of rotating 
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convection in the outermost layer of the convection zone. The 
computations provide a tool for measuring the eddy viscosity 
in the outer solar convection zone. 



2. The Reynolds stress 

We start with the construction of the two components of the 
Reynolds stress tensor, 

Qij = (uj(x, *)«£(*,*)}, (2) 
which describe the angular momentum transport, i.e. 
dQ 

Qr<t> — — v^r-— sin (9 + I ■ Vsm.6, 
or 

dQ 

Qe<j> = ~VT-QQsm9 + I ■ HcosO. (3) 

Due to the terms containing the quantity I, angular momentum 
is transported even in case of rigid rotation, Q = const.. This 
non-diffusive part of the stress, which only exists in case of 
anisotropic turbulence subject to a basic rotation, is known as 
the A-effect (Riidiger 1989). The problem of the angular veloc- 
ity gradient in the outer part of the solar convection zone has 
already been discussed by Gilman & Foukal ( 1 979) and Gailitis 
& Riidiger (1981) with theories based on the A-effect. 

The stress (^) contains three unknown parameters, which 
determine the internal rotation law: the vertical angular mo- 
mentum transport V, the horizontal angular momentum trans- 
port H, and the eddy viscosity v^. If the rotation law is known 
from observations and the A-effect is known from simulations, 
one should be able to measure the eddy viscosity. Once all three 
parameters are known we can compute a theoretical Ward pro- 
file, which can then be compared with the observed one as a 
test. 

In (^|), we have introduced the turbulence intensity, 



(4) 



from the simulations to normalize the A-effect. The functions 
V and H are also known from the simulations, but not the eddy 
viscosity, which we shall choose such that the observations are 
reproduced. With the same normalization for the eddy viscosity 
as for the A-effect, i.e. 



** = V Q> 

the Reynolds stress reads 
I 



(5) 



Qr(p — 

Qe(t> — 



dQ \ 
vr— + VQ sm.0, 
or ) 

dQ 

v-— sin 6 + HQ cos 6 
88 



(6) 



with v as the only free parameter. It is varied between and 1 
in order to reproduce the observed value of 5% for the decrease 
of the rotation rate with radius. It is obvious that large values 
of v will produce negative (positive) cross correlations on the 
northern (southern) hemisphere, while small numbers of v and 
positive values of H are required for the observed (opposite) 



behavior. Indeed, as we have scaled the eddy viscosity with the 
rotation period rather than with the convective turnover time in 
(Q), we must expect a rather small value for v. As the eddy vis- 
cosity in the solar convection zone fixes the value of the Taylor 
number, the Taylor number fixes the amplitude of the merid- 
ional flow, and the meridional flow might fix the solar dynamo, 
it is important to derive the eddy viscosity from observations. 

3. The simulations 

We obtain the functions V and H from the f-plane numeri- 
cal simulations of Chan (2001). These calculations computed 
rotating convection in local pieces of the spherical shell and 
obtained values of the turbulence Reynolds stress at a different 
latitudes. The latitudinal coverage of the cases is dense enough 
for us to obtain analytical fits of the numerical data. The strong 
density stratification is completely included in the model simu- 
lations, as is the energy transport by radiation (gradient of tem- 
perature) and convection (gradient of entropy). 

The resulting A-effect from the Q = 1/2, F/0.25 = 1/8 
case, as tabulated in Chan (2001, Table 3), is given in the Figs. || 
and H showing the data from the simulation vs. the expansions 
of Eqs. (j^) and ([)]). The main feature of the vertical A-effect, 
V, is its negative sign (see Fig. ||) while the horizontal A-effect, 
H, proves to be positive (see Fig. |3|). 

After the theory of the A-effect, the vertical transport can be 
approximated with V oc (vf?) — (u' 7 2 ) for slow rotation, so that 
predominantly vertical turbulence, (u{?) > (u 1 ?), is needed for 
V to assume negative values. As Table 2 in Chan (2001) shows, 
this is indeed the case. 

V (0=1/2) 




0° 

LATITUDE 



Fig. 2. The latitudinal profile of the vertical A-effect, V. Solid 
line: values from Chan (2001), dash-dotted line: The profile 
described by Eq. (M) 



By the signs of both H and V, these results are very promis- 
ing in order to reproduce the negative slope of the rotation rate 
as well as the positive sign of the horizontal cross correlations. 
However, in contrast to expectations, V exhibits a minimum at 
the equator while H is highly concentrated at low latitudes. 

As in Riidiger (1989), the A-effect terms in (§) and (§) are 
written as 



2=0 



H = Y / H (l) sin 2 

1=1 



(7) 
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Fig. 3. The latitudinal profile of the horizontal A-effect, H. 
Solid line: values from Chan (2001), dash-dotted line: The pro- 
file described by Eq. 



so that the result of the simulations can be summarized with the 
series expansion 



V = -0.300 + 0.187 sin 2 9 + 0.0158 sin 4 0- 



Fig. 4. The turbulence-originated coefficients in the non- 
diffusive fluxes of angular momentum in the solar convection 
zone from KR93. Solid line: V {0) . Dashed line: V {1) = if (1) . 
Note that is positive for slow rotation 

rotation rate through the supergranulation layer. With the ex- 
pressions (ph, we solve the Reynolds equation, 



-0.0337 sin b 6», 



and 

H = 0.727 sin 6 9, 
so that 

y(°) = -0.30, 
V {2) = 0.016, 

H (i) = H m = 0i 



V {1) = 0.19, 
V (3) = 0.034, 
ff (3) = 0.73 



(8) 



(9) 



(10) 



P 



■ 5r + («.v)« 



-Vp - V • (pQ) + P g, 



(12) 



results. These results are very interesting insofar as the surface 
effects of the convection zone are included in the model, in 

contrast to the computation of the A-effect for free turbulence f2 — J]] ^n-i(^) 
by Kitchatinov & Riidiger (1993, KR93). In the latter model 
the A-effect is a function of the Coriolis number, 



assuming axisymmetry and under the anelastic approximation, 

V-(pw) = 0. (13) 

In cases where the meridional flow can be neglected, the az- 
imuthal component of Eq. ( |l2| ) is reduced to 

V ■ (rsm9p(u' 4> u')) =0. (14) 

Expanding the rotation rate in terms of Legendre functions, 

P^\cos6) 



sinf 



(15) 



Eq. (J14J) can be solved analytically for a thin layer with stress- 
free boundary conditions, 



n* = 4tt- 



Pn 



(ID 



Sr<t> 



0. 



(16) 



where r corr is the convective turnover time and P rot the rotation 
period. The result is shown in Fig. || For slow rotation, f2* <C 
1, positive values result for and = is small but 
negative. On the other hand, for fast rotation one finds < 
and = > similar to the numbers in Figs. ^ and 

|. In the supergranulation layer of the solar convection zone, 
however, one cannot apply the results for ft* ^> 1. 

The values in (|o|) differ considerably from the predictions 
of the KR93 theory indicating that the latter is invalid for the 
outermost layers of the convection zone, where surface effects 
such as the presence of a boundary, strong density stratification, 
and radiation transport are essential. 

4. The solution 

We are looking for the eddy viscosity parameter v. It must be 
positive and it must reproduce the 5% radial decrease of the 



as described in detail by Riidiger (1989). The boundary condi- 
tions, 



VU) Q = Vq VLl>2 = V2, 



(17) 



where V n now are the components of the function V expanded 
in terms of orthogonal polynomials, lead to 



Vo = V w + 0.8V {1} 
Vi ~ -0.02. 



0.069F (2 



-0.15, 



(18) 



Inserting the slope uj' q ~ —1.4 from the observed rotation law 
in Eqs. tfTfy and @, we obtain 



0.11 



for the eddy viscosity. 



(19) 



The derivation of the result (19) is only valid for a very 



thin surface layer. To treat a layer of finite depth, we solve Eq. 
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(12) numerically with v as an input parameter. At the lower 
boundary, at 2^=0.95, we prescribe the rotation law, 



n(9) = tf (0.7 + 0.3sin 2 



(20) 



to impose the observed latitudinal shear. The upper boundary 
is stress-free. The density stratification is taken from a stan- 
dard solar model (Ahrens et al. 1992). The code used is the 
same as in Kiiker & Stix (2001). We have made runs with and 
without the meridional flow included and found no significant 
difference between the results. The results shown in Fig. || were 
obtained with the meridional flow included. 

Indeed a negative slope of the rotation law results for all 
latitudes. Its amplitude grows with decreasing viscosity. For 



0.1 



(21) 



the results in the Fig. |5] comply with the observations. Again 
the dfi/dr is negative for all latitudes. 

5. Eddy viscosity and magnetic Prandtl number 

As there are no observations about (say) the decay of large- 
scale vortices at the solar surface, we have no direct informa- 
tion about the amplitude of the eddy viscosity. The only way to 
determine it is the study of the internal differential rotation of 
the Sun. We have shown that for v ~ 0. 1 the observed radial 
gradients of the outer solar rotation law can be reproduced. It 
might be interesting to use this value to estimate the eddy vis- 
cosity amplitude in the outermost layers of the solar convection 
zone. With a RMS value of about 200 m/s for the velocity fluc- 
tuations (see Fig. ph, 



iT ~ 1.5 x 10 13 cm 2 /s 
follows from Eq. (|j). 



(22) 




Fig. 6. The turbulence intensity in the solar supergranulation 
layer after Ahrens et al. (1992) 



sunspot decay (~ 10 11 cm 2 /s). The dispersal of large-scale pat- 
terns in the surface magnetic flux, however indicates a consid- 
erably larger value of 6 x 10 12 cm 2 /s for the eddy magnetic 
diffusivity (Sheeley 1992), hence a turbulent magnetic Prandtl 
number of about 2.5. 

Large values of the magnetic Prandtl number are not quite 
unreasonable. It is shown in Riidiger (1989) that the turbulent 
Prandtl number runs with the inverse microscopic Prandtl num- 
ber, which takes the value of 0.01 for the solar plasma and thus 
would indeed lead to values of about 100 for the turbulent mag- 
netic Prandtl number^. 

6. The Ward profile 

An important test is whether such viscosity values would gen- 
erate the positive Ward profile. This is indeed the case. At low 
northern latitudes the horizontal cross correlation is positive, 
but we do not have any information about the higher latitudes. 
From the stress (ph, we find for the Ward profile, 



H cos ( 



(23) 



The results for various values of v are shown in Fig. 0. Note that 
the covariance is indeed positive (in the northern hemisphere) 
at low latitudes, and slightly negative at high latitudes, from 
where we do not have any reliable data. 






In stellar magnetohydrodynamics, the turbulent magnetic 
Prandtl number, ^m=v^lrr, is an important parameter. It is of- 
ten assumed to be of order unity, but this is not finally clear. 
The advection-dominated solar dynamo, e.g., requires an eddy 
magnetic diffusivity smaller than 10 12 cm 2 /s (Choudhuri et 
al. 1995, Dikpati & Charbonneau 1999, Kiiker et al. 2001), 
hence after (^2|) a turbulent magnetic Prandtl number greater 
than 15. The value increases to Pm ~ 150 if - as it is sug- 
gestive - the eddy magnetic diffusivity, rjT, is derived from the 



Fig. 7. The Ward profile W for v = 0.3, 0.1, 0.03 (from top to 
bottom). The correlation is slightly negative in the high north- 
ern latitudes but it is positive as observed in low northern lati- 
tudes 



1 The considerations in Riidiger (1989) only concern the Prandtl 
number rather than the magnetic Prandtl number but the expressions 
are very similar in both cases 
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Fig. 5. Rotation rate at the equator, 15, 30, 45, and 60 degrees latitude (from top to bottom) in the outermost region of the solar 
convection zone for v = 0.3 (left), v — 0.1 (middle), v = 0.03 (right) 



7. The meridional flow 

Doppler measurements find a surface meridional flow that is di- 
rected polewards with maximum speeds between 10 and 15 m/s 
(Komm et al. 1993). Models of solar differential rotation based 
on the KR93 Reynolds stress usually fail to produce this sur- 
face flow, though a corresponding flow cell exists in the bulk 
of the convection zone. It is, however, superseeded by a sec- 
ond flow cell in the surface layers with opposite flow direction, 
hence producing a surface flow towards the equator. This su- 
perficial flow cell is driven by the positive radial gradient of the 
rotation rate, which is due to the short convective turnover time 
(and thus small Coriolis number) in the supergranulation layer. 
Since the meridional flow is mainly driven by the gradient of 
the rotation rate in z-direction, the effect on the flow pattern is 
profound. With the A-effect from the Chan (2001) simulations, 
both the positive rotational shear and the additional flow cell 
vanish. Figure ^] shows the maximum surface flow speed as a 
function of the viscosity parameter v. Contrary to Kitchatinov 
& Riidiger (1999) and Kiiker & Stix (2001), the surface flow 
is directed towards the poles, and (again) for the v = 0.1 case 
(which best reproduces the radial shear) the flow speed lies in 
the observed range. 
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Fig. 8. The amplitude of the poleward directed meridional flow 
at the solar surface 



8. Conclusions 

The theory of turbulent angular momentum transport in stel- 
lar convection zones based on the Second Order Correlation 
Approximation (SOCA) in KR93 predicts a positive vertical 
and vanishing horizontal A-effect in the solar supergranulation 
layer, and negative vertical as well as positive horizontal A- 
effect in the bulk of the solar convection zone. Solutions of the 
Reynolds equation with the stress tensor from KR93 therefore 
reproduce the observed variation of the rotation rate with lati- 
tude remarkably well, but lack the decrease of the rotation rate 
with increasing radius in the outermost part of the convection 
zone. 

Simulations of rotating convection in the upper part of the 
solar convection zone show a strong and positive horizontal 
and a negative vertical A-effect. With the A-coefficients as 
derived from the simulations, the solutions show a negative 
shear, d^l/dr < 0, of the observed amplitude when a value 
of 1.5 x 10 13 cm 2 /s is chosen for the eddy viscosity. As there 
is no way to directly measure the eddy viscosity, this is the only 
method to derive its value from observations. 

As an independent test, we have computed theoretical Ward 
profiles. With the A-effect from KR93, the horizontal cross- 
correlation is always negative because the horizontal A-effect 
vanishes in the surface layer. With the large positive value of H 
from the simulations, on the other hand, Qg^ is always positive 
at low latitudes, as observed, and the amplitudes are in good 
agreement as well. 

In models of the solar differential rotation with positive ra- 
dial A-effect in the outermost layers of the convection zone, the 
radial shear is always positive in that layer, and the surface gas 
flow is directed towards the equator, both in contradiction to 
the observations. The negative radial A-effect derived from the 
Chan (2001) simulations removes both these contradictions by 
maintaining a negative gradient of the rotation rate, which in 
turn drives the surface flow towards the poles. 

We conclude that the KR93 expressions for the Reynolds 
stress are invalid in the outermost part of the solar convec- 
tion zone. Possible reasons are the proximinty of the outer 
boundary, the increasing importance of radiative energy trans- 
port with decreasing depth, and the neglect of the inherent 
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anisotropy of turbulent convection, which is driven by a large- 
scale entropy gradient rather than a random force. 
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